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ABSTRACT 

Galaxies which fall into clusters as part of the same infall halo can retain correla- 
tions due to their shared origin. N-body simulations are used to study properties of 
such galaxy subgroups within clusters, including their richnesses and prevalence. The 
sizes, densities and velocity dispersions of all subgroups with > 8 galaxies are found 
and compared to those of the host clusters. The largest galaxy subgroup provides a 
preferred direction in the cluster and is compared to other preferred directions in the 
cluster. Scatter in cluster mass measurements (via five observables) , along ~ 96 lines 
of sight, is compared to the relation of the line of sight to this preferred direction: 
scatter in cluster velocity di spersion measurements sho w the strongest correlation. 
The Dressler-Shectman test ([Dressier fc Shectman|[l988l ) is an observational method 
to detect cluster substructure. This test is applied to the cluster sample to see whether 
the substructure it identifies is related to these subgroups. The results for any specific 
line of sight seem noisy, however, clusters with large subgroups tend to have a higher 
fraction of lines of sight where the test detects substructure. 
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1 INTRODUCTION 



Galaxy clusters are of interest for many reasons 
for some recent summaries and reviews see e.; 
2005T): 



Voit 



Kravtsov et"aj l|2009l ); iMvers et all ((200' 



Allen. Evrard fc Mantzl ( 201 ll ) and the proceedings from the 
2008 Enrico Fermi School on cluster^]). They are the largest 
virialized objects in the universe and as such extremely sen- 
sitive to the growth rate and mass density of the universe. 
They host unique astrophysical phenomena and the most 
luminous galaxies in the universe and are more generally a 
special environment for transforming galaxies they enclose. 
They were first discovered as clusters of galaxies (hence their 
name); but they are also deep dark matter potential wells 
which gravitationally lens, and are full of hot gas that emits 
in the X-ray and scatters cosmic m icrowave background 
photons (jSunvaev fc Zel'dovicbll 19721 ). All of these different 
properties can be combined to get a fuller understanding of 
the processes within and affecting galaxy clusters, and their 
evolution and observable properties. 

Our interest here is in the galaxies in galaxy clusters, 
specifically those which joined the cluster as part of a shared 
infalling halo. These will be called galaxy subgroups hereon 
and are the focus of this note. (Galaxies falling in from the 
same filament also are expected to have correlations but will 
not be considered here.) Using N-body simulations with re- 



1 http:/ /people. sissa.it/~lapi/efs08_wcbsite/ 
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solved subhalofl identified as galaxies, we study statistical 
properties of these galaxy subgroups for over 200 clusters. 
We measure their quantity, average properties, persistence, 
relation to some intrinsic cluster properties and mass ob- 
servables, and how well they can be detected via one classic 
observational test. 

Some corollaries of the presence of these sub- 



groups are already known dKnebe, Gill & Gibsonl 


2004; 


Knebe et al 2004 Mcintosh et all 


2008: AnEulo et all 


2009; 


bi & Md |200St ISimha et al 200G 


: IWetzel, Cohn & White 


2009: iKlimentowski et all I2010TI. 


For example, satellite- 



satellite mergers within clusters and groups are not uncom- 
mon: galaxies which are group members when they fall in 
often merge with the center of their original infall group. 

Groups of galaxies in clusters which are moving to- 
gether or are more dense in space are sometimes referred 
to as substructur es^ X-ray cluster gas also can exhibit 
substr ucture (e.g. iRichstone. Loeb fc Turner! l|l992h : iBuotd 
(2002) for early work and a review), as can the cluster 
dark matter, e.g. mapped through gravitational lensing (e.g. 
for a study on its eff ects on strong lensing of clusters see 
iHennawi et a 3 l|2007l )). These gas and dark matter sub- 



2 overdensities 



phase space using the algorithm of 



iDiemand. Kuhlen fc Madaul {2006ft . sec below 
3 The term substructure can also refer to subhalos correspond- 
ing to the ga l axies t hemselves in simula tions, early work includes 



iKlypin et all dl999h ; iMoore et all <1999h . 
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structures presumably have some relationship to the sub- 
groups of galaxies of interest here. Similarly, simulations are 
sometimes high enough resolution to resolve subhalos in- 
side of subhalos dWeller Qstriker fc BoddlioO^; IShaw et all 
20071; ISpringel et alll2008l; iGiocoli et alll2009l; iLi fc Mdl2009l ; 
Yang. Mo fc van den Bos ch 2008; Giocoli et al 2010). If one 



identifies these sub-subhalos as galaxies then a relation is im- 
plied between their host subhalos and the galaxy subgroups 
described here. Other related work includes the tracking 
of gr oups of subhal os explicitly inside 8 simulated clus- 
ters (jGill et all[200l l. and the measurements of subgroups 
in Milky Way size simulated hal os, including the ir distri- 
bution and velocity associations l|Li fc Helmill2008h . (Local 
group subgroups often can be identified using more infor- 
mation, such as metallicity or three dimensional positional 
information. There is a large body of literature on the sub- 
ject.) Earlier observational studies of subgroups in clusters in 
particular include analysis of ENACS (ESQ Nearby Cluster 
Surv e y) clu sters and Coma as well (|Gurzadvan fc Mazurd 
1 19981 . 1200 ll ). Previous work is extended here by the joint 
use of subhalos identified as galaxies (which preserves corre- 
lations between galaxies due to their histories), the subhalo 
and halo merging histories, and a large sample of clusters. It 
ex pands upon the subgroup p roperties and examples noted 
in (White. Cohn fc Smitl |2010l ) (WCS hereon), for the same 
simulation. 

As galaxies moving together within a cluster can 
indicate the cluster is not relaxed or, observation- 
ally, that interlopers are projected onto the cluster, 
tests have been designed to detec t them. These tests 
(e.g. the Dressler-Shectman test (|Dressler fc Shectmanl 
Il988h ) have been applie d to in dividual obje ct s (in - 
cluding e.g. iBoschin et all (|2006l ); iGirardi et al l|2008h ; 
Barrena et all d201ll) ) larger cluster surveys (e.j 



Solanes. Salvador-Sole fc Gonzalez-Casadol (I199L.. 

Oegerle fc Hilll (l200ll) ; iBurgett et all ||2004|); iHwang fc Led 
' 20071); iRamella et all f2007h ; iMilvang- Jensen et all | 200; " 



lAguerri fc Sanchez- Janssenl (|2010|); 
and numerical data _fe.g. Icrone 



Einasto et al (201C 
Evrard fc Richstone 



i| 19961 ); ICen] l|l997l) ; iKnebe fc Muellerl l|2000l ); IWCSl) . An 

increase of ga laxy and velocity substructure has been 

corre l ated (lEspino-Briones. Plionis fc Ragone-Figueroal 

l2007l ; iRagone-Figueroa fc Plionis! 120071 ) with more recent 
mergers, higher density environments, and increased clus- 
ter elongation. Pairs of galaxies in clus te rs were studied 
previo usly also, e.g. in Iden Hartogl (|l997l ); [Taylor fc Babull 
(|2004l ). 

The work here is based upon N-body simulations, de- 
scribed in §2. Hydrodynamic simulations, which include gas 
physics and a variety of subgrid prescriptions, are not yet 
available at the volume and resolution considered here. Some 
comparisons of subhalo prop erties with and without hydro- 
dynamics have been made dMaccio et al 2006; Dolag et al 
20081 ; ISaro et alll2008l; I simha et all 120091; Ijiang Jing fc Liii 
20101 ; IKnebe et a3l2010l ; ISchewtschenko fc Maccioll201ll ). In 
some of these comparsions a small fraction of radial orbits 
change (which are themselves a small fraction of subhalo or- 
bits). Direct application to the results here is not straightfor- 
ward as the dark matter subhalo finder used here is different 
than the ones for which the comparisons have been made. 

Simulations and mocks are described in §2, in §3 statis- 
tical properties of the galaxy subgroups are given, §4 gives 



the relation of five mock observational mass measurements 
to properties of galaxy subgroups and to each other, §5 
describes the results of applying the substructure Dressler- 
Shectman test to these clusters, and 86 summarizes. 



2 SIMULATIONS 

A 2048 3 particle, 250h~ 1 Mpc side periodic N-body simu- 
lation box is used, provid ed by Martin White by running 
his TreePM (|Whitd I2OO2I ) code. The cosmological param- 
eters are taken to be h — 0.7, n — 0.95, Q. m = 0.274, 
and erg = 0.8, in accord with a large number of cosmolog- 
ical observations. The simulation has outputs at 45 times 
equally spaced in ln(a) from z = 1 to z = 0. Halo s are 
found using friends-of- friends (FoF) |Davis et al.l ll985) with 
a linking length of b = 0.168 times the mean interparticle 
spacing. Halo masses given below are FoF masses. Our in- 
terest will be in the 243 clusters in the box, i.e. halos with 
M > 1O 14 /i _1 M , at z = 0.1. (This redshift is used be- 
cause some of the mock observations described below rely 
on models which were trained on observational data at this 
redshift.) 

These simulation data were also used in WCS, which 
can be consulted for details of implementations, mock obser- 
vations, and tests beyond those given below. Briefly, galax- 
ies are identified with subhalos and the two words will be 
used interchangeab ly hereon. The subhalos are found v ia the 
FoF6d algorithm of lDiemand. Kuhlen fc Madaul (|2006l ) (the 
specific implementation of their algorithm is described in the 
appendix of WCS). Subhalo i nfall masses can be used to in- 
fer ga laxy luminosities (e.g. IConrov. Wechsler fc Kravtsovl 
(2006)). Here log 10 Mj n f > 11.3 is chosen, which corre- 
sponds to a minimum lumi nosity of ~ 0.2L,, ( at z = 0.1 
this is ~ -18.5 in r band (Blanton et al.l [2003) or stellar 



mass 3 x 10 9 /i _1 M© (|Moster et al.ll2010] )). No luminosity- 
infall mass scatter is inclu ded (some estimates are in e.g. 
Ivan den Bosch et a 3 l|2007l )). In some cases as noted below, 
0.4L, is used as a minimum cut instead. This method of 
luminosity assignments for subhalos gives agreement with 
observations for galaxy clustering, the cluster galaxy lu- 
minosity function, cluster richnesses and the radial cluster 
galaxy profile (see WCS). To find galaxy histories , track- 
ing is as described inlWetzel. Cohn fc Whitel l|2009l ) ; IWCSl ; 
IWetzel fc White] (|2010l ). 

By augmenting the dark matter simulation, WCS mea- 
sured mock observational masses for these clusters via sev- 
eral methods. The relationship between the observational 
masses and the simulated cluster masses, and the form of 
the scatter, are important for understanding the wealth of 
cluster survey data currently in hand and arriving soon 
(see, e.g. the reviews cited above). Five mass measure- 
ment methods from WCS are used here. Two are rich- 
nesses. The first is the iKoester et"al] (|2007l ) MaxBCG al- 
gorithm based upon colorfl The second richness estima- 



4 MaxBCG is one of many algorithms based on a red se- 
quence finder l lGladders fc Yeell20od 120051'). Color assig nments 
are estimated with prescription of Skibba & Shcth (2009), com- 
bined with redshift evolution o f IConrov. Gunn fc Whitel | |2009T) ; 
IConrov. White, fc Gum] fcOlOh ; IConrov, fc Gunnl j2010l V Galax- 
ies are taken to be "red" if they have g — r within 0.05 of the peak 
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tor uses spect roscopy and assi gns cluster membership via 
the criteria of I Yang et all (|2007l )). The third mass measure- 
ment uses Sunyaev-Zel'dovich (SZ) flux or Compton decre- 
ment. Flux is assigned by using halo mass for temperature 
and taking the dark matter particle density proportional 
to the gas density. The flux is then measured within an 
annulus of radius risob, the radius within which the av- 
erage mass is greater than or equal to 180 times back- 
ground density. Tests of thi s appr oximation are given in 
IWhite. Hernauist fc Springell i|2002l ). Weak lensing masses 
are found by using a cluster len s profile of SIS or NFW 
(N avarro. Frenk and Whitejll997t l form and then fitting for 
a velocity dispersion and then mass. The velocity disper- 
sion masses rely upon phase space information to reject 
outliers and include a measured harm onic radius in the 
mass c alculation, based on methods of Ivan Haarlem et al.l 
ill997h;lden Hartog fc Katgert] (ll996l ) ;lKatgert et all (|l996h ; 
iBiviano et al l (|2006l ); IWoitak et all (|2007l . I2009D F1 We use 
the WCS measurements corresponding to radius risob when 
a radius needs to be specified. WCS measured each indi- 
vidual cluster's mass along 96 different lines of sight. This 
ensemble of cluster mass measurements will be used here as 
well, and as in that analysis, lines of sight where a more 
massive cluster has its center within ri80i> are removed. In 
addition, for the work here, lines of sight where either rich- 
ness was < 2 were also discarded. 



3 PERSISTENCE OF GALAXY SUBGROUPS 

The inhomogeneities in cluster galaxy distributions are in 
part historical, as larger halos are built up from the infall 
of smaller halos. When halos fall into larger halos, they be- 
come subhalos themselves, with their central and satellite 
galaxies all now becoming satellites within the new larger 
halo. To study these galaxy subgroups, for every z — 0.1 
cluster galaxjO its infall group, infall time, infall group rich- 
ness (above the luminosity cut) and infall group mass are 
identified. Each infall halo containing more than one galaxy 
results in a separate galaxy subgroup. An example of a 9 
galaxy subgroup is shown in Fig. [1] along with the tracks 
of these galaxies since their host halo fell into the cluster 
6.3 Gyr earlier (i.e. all the tracks are tracks of the subgroup 
members within the cluster itself). At infall these galaxies 
were m a 2.5 x lO 13 ft _1 M halo with 12 galaxies. Triangles 
mark the final galaxy positions at z = 0.1. It can be seen 
that the galaxies have been staying together throughout. 
The dashed circle is risot of the host cluster at z — 0.1, cen- 
tered at the cluster center. The host cluster has 313 galaxies 
and mass 9.2 x 10 14 /i' 1 M Q . 

A large fraction of galaxies within clusters are in these 
subgroups. At our redshift of interest, z ~ 0.1, >40% of 



of the red galaxy g — r distribution specified bv lSkibba fc Shethl 
(2009), for their observed M r , again see WCS for more detail. 
5 We will not use the velocity dispersions based upon 3 a clipping, 
presented in WCS, as these were less well correlated with the FoF 
cluster mass. 

B Only true cluster galaxies are included in this section, interlop- 
ers are considered in §5 and §6 for observational mass estimates 
and substructure estimates. 
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Figure 1. Tracks of subgroup galaxies after infall into a cluster. 
This subgroup fell into the cluster 6.3 Gyr ago, and contained 
12 galaxies at that time. Only tracks of the galaxies surviving to 
z = 0.1 are shown; solid triangles mark their z = 0.1 positions. 
The dashed circle is r 180 i, for the full 9.2 xl0 14 h~ 1 Mo cluster, 
which has 313 galaxies. The center position (0,0,0) is the cluster 
center at Z = 0.1. 



the cluster galax ies ab ove 0.2L* have at least one associ- 
ated companion l|WCSl ). (As mentioned earlier, these sub- 
group galaxies presumably corres pond in some way to the 
sub-s u bhalos mentione d earlier dWeller. Qstriker fc Bode 
| 200, ' 
pOO! 

iGiocoli et all [2010|).) Not all the galaxies which are not 
in subgroups fell in alone: of the galaxies which have 
no companion in the clusters at z = 0.1, ~ 1/5 fell 
in with companions which have disappeared since infall. 
Some of the subgroup galaxies (~ 30% of cluster galax- 
ies) were preprocessed in groups, i.e. fell in from halos with 
Mhalo > 10 13 Mq. The fraction of prepr ocessed clus- 
ter g a laxies is similar to e arlier measurements (jBerrier et all 
2006; iMcGee et all 120091 (the simulated galaxy samples 
are not identical, so exact agreement is not expected). 
As in those papers, a larger fraction of the more massive 
galaxies were processed in group environments. The impor- 
tance of the group environment for galaxy formation has 
been stressed in these papers and in e.g. Zabludoff et a] 



IShaw et all 120071 ; ISpringel et all |200S|; IGiocoli et aJ 
iLi fc Mol 120091 ; lYang. Mo fc van den Boschl 120081 ; 



.Zabludofj (120071), see al s o IZabludoff fc Mulchaev 
200d ): lLi. Yee fc Ellingsonl (|2009h . 



The distribution of subgroup richness per cluster (i.e. 
average fraction of subgroups with given richnesses) is shown 
in Fig. [2] for three different mass regimes (44 clusters in 
each). It appears to have a self-similar distribution, al- 
though there is a sharper cutoff at high richness for the 
lower mass clusters. (This is at least in part likely to be 
a resolution effect, as the maximum number of galaxies in 
a subgroup is limited by the total number of galaxies in 
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Figure 2. Fraction of subgroups, per cluster, with given richness, 
averaged over 44 clusters in three mass ranges. The subgroup 
richnesses, similar to the infall halo masses for clusters, seem to 
have a roughly self-similar distribution, although there is a cutoff 
at high subgroup richness (i.e. > 20) for lower mass clusters. 



the cluster, which is smaller for the lower mass clusters.) 
There is likely some way to connect this subgroup distribu- 
ti on to the infall halo mass distribution, s uch as measured 
De Lucia et all (12004); iGao et al.l (l2004|);lTavlor fc Babull 



as the in- 
masses (see 



|2005l ); iGiocoli et all (|2009l ); iGao et all (|2011l 
fall group richnesses are related in infall 
also IShaw et~a 1 (2007) which compares instantaneous sub 
subhalo masses). 

The richness distribution of the largest subgroup in each 
cluster is shown in Fig. [3] Although the largest galaxy sub- 
group for many clusters has ~ 10 galaxies, much larger 
subgroups do occur. The largest subgroup in our sample 
has 86 galaxies, and is in the second largest cluster (M = 
1.1 x 10 15 /i _1 M Q ). This subgroup had 107 galaxies and host 
halo mass 2.4 x 10 14 /i _1 M© upon infall 2.4 Gyr ago. Divid- 
ing largest subgroup richness by host cluster richness, one 
sees that the largest galaxy subgroup has on average about 
20% of the richness of its host cluster, but that there is a 
wide scatter (Fig. [3] bottom) . The ratio of subgroup rich- 
ness to cluster richness might be expected to be relevant 
because more massive halos tend to form from mergers of 
more massive progenitors, and richness is roughly propor- 
tional to mass. (For reference, a cluster of mass 1O 14 /i _:l M 
has richness of ~ 40 on average.) The particular richness 
fraction of the largest galaxy subgroup is not only deter- 
mined by the infall mass of the cluster progenitor, but also 
the time since infall (i.e. the rate of disappearance of the sub- 
group galaxies and the growth of the cluster richness since 
infall). Similarly, the time since a recent major merger is 
correlated with the size of the largest halo subgroup. (Two 
standard criteria for a major merger are that the smaller 
to larger mass is 1:3 or larger or 1:10 or larger, both were 
considered here.) That is, the larger the subgroup the more 
recent a major merger. The strongest correlation is with the 
fraction of the cluster richness due to the halo subgroup and 
the time of most recent 1:3 merger. (A related correlation 
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Figure 3. Top: for each of the 243 dark matter halos with 
M > 10 14 h~ 1 M Q , the largest number of galaxies remaining from 
the same original infall halo (subgroup richness). Bottom: the 
richness of the largest galaxy subgroup as fraction of its host 
cluster richness. More massive halos are are more likely to have 
more massive subgroups of galaxies. The average value of largest 
subgroup richness over host cluster richness is 0.17. 



was f o und by Espino-Briones. Plionis &: Ragone-Figueroal 



(|2007n : lRagone-Figueroa fc Plionisl l|2007h : thev found that 
measured cluster substructure increased for clusters with re- 
cent mergers.) The median subgroup richness fraction (of 
the host cluster's richness) is about 1/3 for clusters with a 
1:3 merger in the last time step (as required by the major 
merger definition), while for clusters with a 1:3 merger ~ 10 
Gyr ago, the richness fraction is closer to 1/10. This is not 
surprising as after the major merger, the cluster will grow 
in richness, while the subgroup will decrease. 

Instead of the largest subgroup per cluster, one can con- 
sider all "large" (defined to have > 8 galaxies, chosen for 
convenienc^ll) galaxy subgroups. Their frequency per clus- 
ter is shown at the top of Fig. U 118/243 clusters have no 
large subgroup, but 38 clusters have more than one large 
subgroup, and three clusters have 5 large subgroups. Almost 
all (42 of the 44) clusters above 2.5 x 1O 14 /i _1 M contain 
at least one large subgroup. For lower mass host halos, large 
subgroups are more rare (e.g. occurring only in 1/4 of the 
725 halos with mass > 5 x 10 13 /i -1 M Q ), when they occur 
they are a larger fraction of the halo richness. Shown also 
in Fig. [4] bottom, is the fraction of cluster galaxy richness 
in large subgroups: the solid line is richness fraction in large 
subgroups for the ~ 50% of clusters which contain at least 
one large subgroup (peaked around 15-20% but with a long 



7 At z = 0.1, richness 8 corresponds roughly to a halo with M ~ 
2 x lO 13 h _1 M0, but of course the subgroups are remnants of 



larger galaxy subgroups at infall. 
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Figure 4. Top: the number of large (> 8 galaxies) subgroups 
in all clusters. Bottom: (solid line) for the 50% of clusters (M > 
10 14 /i — 1 Mq) with subgroups of > 8 galaxies, the fraction of total 
cluster richness due to these large subgroups, and (dashed line), 
for all clusters, the fraction of total cluster richness due to galaxies 
with no associated galaxies in the cluster, per cluster. 



Figure 5. Average fraction of galaxies remaining for infall groups 
initially with (lines top to bottom at left): >20 galaxies (dot- 
dashed), 15-20 galaxies (long-short dashed) 11-14 galaxies (solid), 
7-10 galaxies (dotted), and 3 galaxies (dashed), in all halos with 
M > 1O 14 /i _1 M0 at z = 0.1. Groups where all galaxies have 
disappeared are not shown. Note the minimum possible fraction 
of galaxies remaining changes with infall richness. 



tail to large fractions). For comparison, the richness fraction 
in all clusters due to all galaxies not in subgroups is shown 
as well. The latter is 60% on average, as mentioned above, 
but has a wide spread. 

These z = 0.1 subgroups are remnants of larger sub- 
groups of galaxies upon infall. Some subgroups fell in long 
ago: over half of the 189 large subgroups fell in more than 
2 Gyrs earlier. The subgroups which have fallen in more re- 
cently tend to have a larger surviving fraction of galaxies. 
For example, the surviving fraction of galaxies per current 
large subgroup has a median value of ~94% for the 36 (cur- 
rently) large subgroups whose halos fell in 0.6 Gyr ago, drop- 
ping to ~ 50% for those 26 falling in 4.9 or more Gyr ago. 
Considering instead the initial richness of infalling groups, 
the remaining fraction of galaxies (if nonzero) is shown for 
subgroups which have 3, 7-10, 11-14, 15-20, and > 20 infall 
galaxies, as a function of infall time, in Fig. [5] Within the 
noise the change in galaxy fraction looks similar for different 
initial galaxy richnesses (for subgroups starting with only 3 
galaxies, the fraction is bounded below by 1/3). Because 
of this galaxy attrition within subgroups, the satellite with 
the largest infall mass (i.e. formerly the central galaxy of 
the largest infalling halo) is not always in the largest sub- 
group. A smaller halo might have a larger remaining galaxy 
fraction. For 70% of the simulation clusters, the largest sub- 
group does have the satellite with the largest M[ n f. Another 
5% of the clusters have the satellite with largest infall mass 
in a large subgroup (i.e. > 8 galaxies) but not the largest 
sub group. When using subhalo abundan ce matching (see, 
e.g. IConrov. Wechsler fc Kravtsovl |2006l )), the galaxy with 



the largest infall mass will be the brightest satellite (up to 
scatter). However, given the above, even in the case of no- 
scatter (i.e. oversimplified) halo abundance matching, the 
brightest satellite will not be a useful way to identify the 
largest subgroup. 

The spatial and velocity distributions of large subgroups 
differ on average from those of their host clusters. As ex- 
pected from their lower richnesses and their membership in 
a larger cluster, subgroup sizes are generally smaller: the 
average distance of subgroup members from the subgroup 
center is about half that of all cluster galaxies from the host 
cluster's center, albeit with scatter|f] (However, one of the 
189 large subgroups has an average radius larger than that 
of its host cluster. The subgroup center is taken to be the 
average position, as not all subgroups have the infall central 
galaxy remaining.) The average distance of subgroup galax- 
ies from their center tends to be larger for older subgroups. 
One possibility is that galaxies close to the subgroup cen- 
ter have more time to merge with the subgroup center and 
disappear, while galaxies further away have more time to be 
pulled away from the center by the host cluster tidal fields. 

In terms of density, modeling the cluster and the sub- 

8 If the subgroup radii are rescaled by M.^^ 3 of their infall ha- 
los, and similarly for the full cluster and its mass, the ratio of 
subgroup to cluster radii becomes very broadly centered around 
unity, that is, further suggesting that the smaller size of the galaxy 
subgroups is partially due to their smaller size at infall. I thank 
M. George for suggesting this and the rescaling below. 
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Figure 6. Left: The "density" (number of galaxies over ellipsoidal 
volume) of large (N g > 8) subgroups divided by its counterpart 
for all galaxies in their host clusters. Right: velocity dispersion of 
galaxies in large subgroups relative that of all galaxies in their 
host clusters. For both left and right, there can be more than one 
large subgroup in the same host cluster, as shown in the top panel 
of Fig.Q] 



groups as ellipsoids and then calculating density as the num- 
ber of galaxies divided by respective ellipsoidal volumes, 
the large galaxy subgroups have a median (average) den- 
sity 30%(80%) times larger than clusters they live in. The 
broad distribution is shown at left in Fig. [6]-almost 1/3 of 
the large subgroups have densities which are more than dou- 
ble that of their host cluster. In terms of shape, the median 
value of the (minimum axis/maximum axis) for large sub- 
groups divided by that for their host cluster is slightly below 
1 (~ 4/5), but there is a tail out to values> 3. That is, the 
majority of subgroups (~ 60%) tend to be less round than 
their host clusters. Subgroups which fell in earlier tend to 
be more elongated compared to those which fell in later (the 
average short/long axis ratio for subgroups falling in at the 
last time step is ~0.4, while for subgroups which fell in over 
6 Gyr ago, the average short/long axis ratio is ~ 1/4). It 
might be expected that this elongation is due to tidal forces 
within the larger cluster stretching out the subgroup after 
its infall. However, there isn't a strong signal for alignment 
of the long axis of the large subgroups with their direction 
of motion within the cluster. 

In terms of velocities, the subgroups tend to have 
smaller (~70%) velocity dispersions (relative to the average 
subgroup velocity) compared to their host cluster's coun- 
terpart, i.e. they are slightly colder(f] The ratio of velocity 
dispersions of large subgroups to those of their host cluster 
is shown in Fig. HJ] The subgroups which fell in earlier tend 
to have larger velocity dispersions, both in raw numbers and 
(less strongly) relative to that of their host halo. 

The radial distribution of the centers of the largest sub- 
group per cluster (rescaled by the rms distance of the cluster 
galaxies from the cluster center) lies alon g an N FW profile, 
just as the satellites in the simulation do (|WCS() and as has 
been seen observationally for clust er satellites in stacked pro- 
files l|Lin. Mohr fc Stanfordll2004l ). It seems the outskirts of 



9 Scaling the subgroup velocities by their infall halo M. 1 ^ 3 and 
similarly for the cluster hosts actually gives a median velocity 
dispersion which is higher for the subgroups, that is, the sub- 
groups are colder than their hosts but hotter than isolated halos 
of similar mass. 
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Figure 7. Left: The distribution of inner products of average po- 
sition direction of the largest subgroup (top) and all large (> 8 
galaxies) subgroups (bottom) of galaxies, with the long axis di- 
rection of the M > 10 14 ft — 1 Mq clusters, as determined by the 
moment of inertia of its dark matter particles in its FoF halo. 
Right: The distribution of inner products of the average velocity 
direction of galaxies within the largest subgroup (top) within the 
cluster and all large (> 8 galaxies) subgroups (bottom) of galax- 
ies, with the long axis direction of the M > 10 14 h~ 1 Mq clusters. 
The large (> 8 galaxies) subgroups and largest subgroup per clus- 
ter tend to both lie along the long axis of the cluster and to be 
moving, in average, along this axis. 



the stacked subgroup position profile are fit better by a lower 
concentration (i.e. 3), but the statistics are very noisy. 

Clusters have several preferred directions, which often 
are related: halos and thus their galaxies tend to fall 
in along filaments. Filaments connecting clusters (i.e. 
usually the largest filament in the cluster) tend to lie 
along the cluster major axis, and the central halo is often 
aligned with one or the other or both (references for 
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ISales et all l|2007l) ; iFaltenbacher et all i|200Sf ); |Paz et all 
l|201ll )), an d has been characteriz ed anal yti cally as well 
(see, e.g. IPennarrubia fc Benson! l|2005h : iLee fc Kang I 
(2006)). As the largest subgroups tend to originate from 
massive halos, they might be expected to lie along the long 
axis of the cluster, i.e. the cluster's longest ellipsoidal axis. 
This is borne out in the simulated clusters: large subgroups 
and the largest subgroup per cluster both tend to lie along 
the long axis of the cluster and to be moving (i.e. have 
average velocity) along this axis (Fig. [7J. The clusters 
with longer large axes tend to have larger subgroups, and 
subgroups further away from the center, however both of 
these may be due to correlations between cluster mass and 
subgroup richness. 

Concentration is often strongly c orrelated with clus- 
ter his t ory and other properties (e.g. see lJeeson-Daniel et all 
l|201lft ; ISkibba fc Macciof (|201ll )). In our sample, the cluster 
concentrations and deviations from the average concentra- 
tion for a given mass were correlated with several quanti- 
ties (the size of the largest subgroup, the deviation of the 
size of the largest subgroup from the average for that mass, 
the fractional richness in large subgroups, and the deviation 
of the fractional subgroup richness in large subgroups from 
the average). However, these correlations were weak, ~ 10%- 
20%. They were in the expected sense, that larger subgroups 
tended to be in clusters with smaller concentrations. 

To summarize, the majority of galaxies in clusters have 
no accompanying galaxies remaining from their infall host 
halo, but galaxies in galaxy subgroups still contribute a sig- 
nificant fraction of cluster richness. Galaxies in large sub- 
groups (> 8 galaxies) tend to have smaller relative velocities 
and higher relative densities than their hosts, even if their 
infall was several Gyr earlier. 



4 SUBGROUPS AND MASSES 

Observationally, a cluster's galaxy population is often used 
to find the cluster's mass, for instance via richness and veloc- 
ity dispersions. The large subgroups described in the previ- 
ous section are deviations from an idealized cluster of virial- 
ized galaxies within an isolated spherical dark matter poten- 
tial. As such, they might affect the use of galaxies as tracers 
of cluster mass. In this section the correlation between scat- 
ter in cluster mass measurements is compared to properties 
of the large galaxy subgroups, when present. 

In WCS, the dark matter simulations were augmented 
as described in §2, and mass measurements were made 
along 96 different lines of sight for every cluster, using 
6 different techniques, for massive clusters (Afigot > 2 x 
10 14 /i _1 Mq). In that paper, some correlation between mass 
scatter and line of sight substructure (via the Dressler- 
Shechtman (|Dressler fc Shectman|[l988l ) test, see below) or 
large or small halos contributing gal axies outside the clu ster 
was found. In the same simulations iNoh fc Cohnl (|201lt ) of- 
ten found an increase in observable mass scatter when view- 
ing along the plane containing the most (filamentary) mass 
in a 10 h~ 1 Mpc radius sphere around the cluster, or along 
the long axis of the cluster. Substructure and observational 
mass scatter have bee n considered in othe r contexts, e.g., in 
X-ray measurements, Ijeltema et all l|200ct ) found mass scat- 
ter correlations with X-ray substructure. 



It should be noted that since these mass measurements 
are only within the 250 h~ 1 Mpc box, scatter due to the 
full line of sight is underestimated. Both weak lensing and 
Compton decrement measurements can easily have contri- 
butions to scatter from well outside these scales. This addi- 
tional uncorrclated scatter would decrease the correlation of 
the observed scatter with intrinsic cluster quantities studied 
here. 

Taking the largest subgroup of associated galaxies for 
each cluster, we calculated the projection of its position on 
viewing direction for ~96 directions 10 !. For many of the clus- 
ters, observed masses increased when measurements of mass 
were along the largest subgroup's axis, i.e. looking along the 
direction where the subgroup is almost directly behind or in 
front of the cluster center. 

One example in shown in Fig. [5] This shows the frac- 
tional mass scatter for a cluster of mass 1.5 x 10 14 /i _1 M© 
and 66 galaxies, for five mass measurements. Its largest sub- 
group fell in 0.61 Gyr ago and has 12 galaxies. Most of the 
observational mass estimates for this cluster increase as the 
line of sight tends to the axis collinear with the largest sub- 
group (note this is also related to the long axis of the cluster, 
as in Fig. [7] and tends to be in the plane containing most 
of the filamentary mass feeding the cluster) . Correlation co- 
efficients are calculated using all points. The correlation of 
scatter in mass is largest for the two measures of richness 
and velocity dispersions, as might be expected. (Again, the 
large SZ correlation is likely overestimated, as is the corre- 
lation for weak lensing, because the small box size does not 
include uncorrelated scatter from larger scales where lens- 
ing and SZ are also sensitive). Also shown at lower right is a 
map of the cluster galaxies with the largest subgroup galaxy 
positions noted (centered 2.2 h~ x Mpc from the central clus- 
ter galaxy). Several of the clusters with recent mergers still 
have the subgroups on the "edge" of the cluster, however, 
the dark matter density in the region between these galaxies 
and the rest of the cluster galaxies is all above the threshold 
set by the b — 0.168 linking length (roughly 100/%). 

The range of line of sight mass scatter varies by clus- 
ter and mass measurement. On a cluster by cluster basis, 
the 65th-35th percentiles of mass scatter (i.e. the 65th - 
35th percentiles of (Af mcas — Mt fU e) /Af truc , roughly the range 
of middle 1/3 of mass scatter) on average is ~ 20%, 10%, 
5%, 30% and 15% for MaxBCG richness, phase richness, SZ 
mass, velocity dispersion mass and weak lensing mass re- 
spectivelvF 1 ! The maximum to minimum mass scatter per 
cluster due to line of sight is on average about 10 times 
this (for SZ it is ~ x30 ). The width of the 65%-35% mass 
scatter per cluster increases with fraction of richness in the 
cluster's largest subgroup for SZ and phase richness, and de- 
creases for weak lensing. The maximum to minimum widths 
are correlated with the fraction of cluster richness in the 



10 Again, lines of sight with a more massive cluster within riggb 
or with either richness < 2 are discarded. 

11 This scatter is due only to changing the cluster line of sight. 
For all clusters, rather than one cluster along different lines of 
sight, this width approximately doubles, except that for velocity 
dispersions, which only increases slightly. As SZ scatte r is often 
closer to 20% in larger boxes, e.g.. lCohn fc White (20Q3), its local 
scatter measured here is expected to double when uncorrelated 
structures are included. 



© 2011 RAS, MNRAS OOO.miHl 



8 J.D. Cohn 




-0.2 0.2 0.4 -16 -15 -14 -1 -0.5 0.5 1 

4M / M i™ Mpc/h Vel 

cluster 112. M= 1.6 xlO" h"' M„ 



Figure 8. Correlation of mass scatter with projection of dis- 
tance to largest subgroup along line of sight, for a cluster with 
mass 1.5 X 10 14 /i — 1 Mq and 66 galaxies, for (left to right, top to 
bottom) red galaxy richness, richness via phase space, SZ, weak 
lensing and velocity dispersions. This cluster has 66 galaxies, the 
largest subgroup at z = 0.1 fell in 0.61 Gyr ago and has 12 galaxies 
and the second largest subgroup has 6 galaxies and fell in almost 
5 Gyr ago. The center of the largest subgroup is 2.2 h~ 1 Mpc 
from the cluster central galaxy. Correlations shown in each panel 
are calculated including all points shown. At lower right, pro- 
jected onto one plane, are the cluster galaxies (pinwheels), with 
the largest subgroup galaxies represented by large filled in points. 



largest subgroup when the subgroups are at least 20% of 
the cluster's richness (69/243 of the clusters). 

For these five different observable masses, one can again 
measure the distribution of correlations between observa- 
tional mass scatter and projection on the direction of the 
largest subgroup of the line of sight. These correlation co- 
efficients are shown in Fig. [9] Although not all clusters had 
correlations between the line of sight to the largest subgroup 
and mass scatter, several did. Richness and velocity disper- 
sions tend to be enhanced when the line of sight tended to 
the direction of the largest subgroup in the cluster, as in the 
example in Fig. [8] These trends can be compared to the cor- 
relations of mass scatter with the long axis of the cluster and 
with measurements through or within the filamentary plan e 
around the cluster, studied in IWCSl ; iNoh fc Cohnl l|201lf ). 
Correlations of mass scatter with respect to the direction of 
the largest subgroup are comparable but smaller for the two 
richness measurements and much smaller for SZ and weak 
lensing. However, the average correlation between mass scat- 
ter and direction of observation relative to the largest sub- 
group is much larger for velocity dispersions than the corre- 
lations found using other cluster preferred axes. 

Combining all the clusters, the correlation between 
mass scatter and direction to largest subgroup, or largest 
subgroup if one is present with > 8 galaxies, is much weaker. 



Figure 9. For each of the 243 dark matter halos with M > 
W 14 h~ 1 Mq, the correlation between the distance to the largest 
associated subgroup of galaxies, projected on the line of sight, 
and the mass scatter. Measuring mass along the direction of the 
largest subgroup of galaxies often produces increased scatter in 
masses based upon richness and velocity dispersions. Smaller cor- 
relations are seen for SZ and weak lensing (it should be noted 
that for these two measurements that the simulation box size 
is relatively small compared to the line of sight distances which 
contribute to their scatters). 

This is presumably due to the tails into negative correlations 
which are seen for some clusters in Fig. [9] One possible cause 
for weak or opposite sign correlations between mass scatter 
and the line of sight for measurement versus direction to 
largest subgroup is the presence of more than one large sub- 
group in many clusters (e.g. Fig. [4]). 

Another quantity was explored as well, the correlation 
of scatter with projection onto the position of the satellite 
galaxy with largest Mj n f. This galaxy would be, in the ab- 
sence of luminosity-infall mass scatter, the brightest non- 
central galaxy. The relation between the direction to this 
galaxy and the line of sight of mass measurement was not 
as strongly correlated with the mass scatter (i.e. the distri- 
bution of correlation coefficients for all the clusters tended 
to be centered on zero). 

To summarize, mass scatters in cluster measurements, 
for several different methods, are correlated with relative 
directions between the line of sight and the axis pointing 

12 Many of the mass scatters are correlated with each other, due 
to their origin in pro perties in and around the cluster, a n d hav e 
been studied in e.g. Word et al.l ||2008|); ICohn fc White) J200ST) ; 
iMeneghetti et al.l ll2010h:IStanek et al.l teOirjMWCSl. As d i scusse d 
in iNord et alj feOOSl ); iRvkoff et all feOOST ); Istanek et alj fculd) ; 
IWCSl correlated mass scatter can not only cause underestimation 
of mass errors for single objects, but can also introduce biases 
when measurements are stacked. 
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towards the largest cluster galaxy subgroup, similar to cor- 
relations found earlier between mass scatter and directions 
of observation relative to the long axis of the cluster and 
to the filamentary plane around the cluster. Velocity disper- 
sion mass scatter is more correlated with the direction to the 
largest galaxy subgroup than with the long axis or filament 
plane directions. 



5 FINDING SUBGROUPS 
OBSERVATIONALLY 

As seen in earlier sections, large galaxy subgroups are 
frequent in massive clusters, and the relation of their 
position relative to the line of sight often correlates with 
mass scatter, particularly velocity dispersions. The ability 
to detect these subgroups of associated galaxies thus may 
be useful both for their own interest and for studies aimed 
at galaxy cluster masses. The presence or absence of galaxy 
substructures (galaxies close in space with smaller relative 
velocities) is also sometimes used to estimate whether a 
cluster is relaxed or not, for example to allow the assump- 
tion of hydrostatic equilibrium for X-ray analyses. Several 
techniques e xist for finding galax y subs t ructure within 
cluste rs fe.g. | Dressler fc Shectmanl (Il988h; iPinknev et alJ 



1996); iKnebe fc Muellerl (|2000f ); iGurzadvan fc Mazurd 
20011 ); iHou et al.l (|2009h and those noted earlier). These 
have been used in several observational and numerical 



studies 



m 

inc luding 



Crone. Evrard fc Richstonel dl996l) 



ill999h; IKnebe fe Muellerl d200Ch: lOegerle fc Hill (l200ll) 



the local counterparts include only the N nn nearest galax- 
ies to the chosen galaxy, in the plane of the sky. The sum 
A = X] S ai $9 al over au cluster galaxies is used to estimate 
whether substructure is present. In the original formulation, 
A is compared with random velocity shufflings in the same 
cluster. A low fraction of shuffled directions giving A larger 
than the unshu ffled A (P < . 05) is considered a substruc- 
ture detection. IPinknev et al.l (| 19961 ) take instead N nn as 
the square root o f the t otal number of cluster galaxies and 
IKnebe fc Muellerl ((2000) take N nn = 25. The latter found 
that A/iVg a i > 1.4 was also a useful criterion to detect sub- 
structure, using dark matter particles, other studies have 
found that a mass dependent thres hold in A is more accurate 
ijRagone-Figueroa fc Plionisll2007l ). In addition to choosing 
the number of nearest neighbors, N nn , it is also necessary to 
specify a luminosity threshold when applying the test, and 
to select which galaxies are in the cluster. 

For application to our mock catalogues, a priori knowl- 
edge of true cluster galaxies cannot be assumed if observa- 
tions are to be mimicked. Cluster membership for galaxies 
must instead be determined using observational techniques. 
The original DS test used 3 — a clipping, but more sophis- 
ticated methods exist for id e ntifyi n g cluster members (for 
examp l e. Ivan Haarlem et all (1 19971); den Hartog fc KatgertJ 
||1996I) ; iBiviano et all (|200d ) Hwoitak et all |2007l , \200m 
which take into account both the line of sight velocities and 
the distance from the estimated cluster center. We use the 
variant of these methods as outlined in WCS to identify 
cluster galaxies. We considered Mi n f cuts corresponding to 
ICenl ^ (|l997l); [jglanes. Salvador-Sole fc Gonzalez-Casadd L > 0.2, 0.4L«, and N nn = ^fW g , 10, 25, iV 9 /6. This last was 

aimed at getting at the natural subgroup size implied by 
Fig.El 

The presence of cluster substructure seemed to depend 
strongly upon viewing angle. Although for some clusters 
substructure was detected more often when the line of sight 
was perpendicular to the largest subgroup, for other clus- 
ters no correlation was apparent between the direction to 
the largest subgroup and the amount of substructure found. 
Taking in principle the cleanest cases, the 87 clusters with 
only one 1 large subgroup, substructure tended to be more 
often detected if the subgroup was perpendicular to the line 
of sight. In this case the subgroup was on the 'side' of the 
cluster in the plane of the sky. But even here, ~ 1/4 of the 
clusters had substructure detection increasing as the axis to 
the substructure became more closely aligned to the line of 
sight. 

Rather than considering substructure detection along 
any particular line of sight, one can instead consider the 
number of lines of sight with substructure, as a property 
of a given cluster. Fig. [TU] shows the distribution, for all the 
clusters, of the fraction of lines of sight where substructure is 
found, by both measures, using N nn — \/N g , and minimum 
infall mass M- ln i > 11.69 (to roughly correspond to L > 
0.4L*). The dotted line corresponds to the fraction of lines of 
sight which meet the substructure detection requirement of 
^sub < 0.05. The numbers of lines of sight with substructure 
with the two measures tended to be correlated with each 
other. Hardly any clusters had no substructure along any 
line of sight. Many clusters had 10% of their lines of sight 
showing substructure, and a handful had over half of their 
lines of sight showing substructure. 

The results depended upon the parameters and sub- 



iBurgett et all (120041); iBoschin et all feOOrj); iHwang fc Led 



(120071); iRamella et all ll2007h; Idrardi et all [|2008h; 
iMilvang- Jensen et al.l (|200ST> ; lAguerri fc Sanchez- Janssenl 
(|2010l ); lEinasto et all (|2010l) ; lwCSl ; iBarrena et all (l201ll). 

Here we apply the classic test by Dressier fc Shectmanl 

j|l988l FI to see how its substructure detection criterion re- 
lates to the presence of large subgroups. The study here 
differs from previous work by focussing on a specific sort of 
substructure present in the clusters, substructure in galaxies 
coming from the same infall group. It also differs from some 
previous works as it uses subhalos rather than dark matter 
particles as galaxies in applying the test, and thus can ap- 
ply interloper methods to identify cluster galaxies in closer 
analogy to observations. (That is, just as in WCS but not 
in most previous work, whole subhalos, with dark matter 
central position and mean velocity, are used as galaxies for 
the tests, rather than random dark matter particles.) 

The original DS test checked for collections of A r nn + 1 = 
11 galaxies whose relative velocities are less likely than ran- 
dom within the cluster ( "random" is calculated by shuffling 
the velocities of the galaxies in the cluster). One calculates, 
centering on each galaxy, 



c2 Nnn r/_ _\ 

°gal = IT ("local - V) 

O 



+ (O local — CT ) 



(1) 



In the rest frame of the cluster, v and a are the mean veloc- 
ity and velocity dispersion of all the cluster galaxies, while 



13 IPinknev et al. (1996) find that DS is the most sensitive to 
substructure, but that it also has the highest detection of non- 
substructures due to elongation and velocity dispersion gradients. 
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considered, then substructure was more often detected for 
clusters whose largest subgroup was farthest from the clus- 
ter center in three dimensions (with a stronger correlation 
with increasing N nn in the substructure test). 

In summary, the Dressler-Shectman test did tend to 
detect substructure more often in clusters with larger sub- 
groups, but the results varied strongly between lines of sight. 
The increased frequency of sightlines which detect subgroups 
as the subgroup size increases is encouraging. However, the 
lack of substructure as defined by this test does not neces- 
sarily mean the substructure is absent. This suggests using 
caution when taking the results of the DS test to measure 
whether substructure is present or whether the cluster is in 
hydrodynamic equilibrium (e.g., for the application of X-ray 
mass estimates). 



Figure 10. Fractions of lines of sight, per cluster, which have 
substructure by the Dressler-Shectman test (described in text). 
The number of galaxies per test subgroup is taken to be the square 
root of the number of cluster galaxies (found via method described 
in text as well) with L > 0.4L*. The solid line corresponds to the 
number of clusters where substructure is detected along a given 
fraction of sightlines, via the condition A/JV ga i > 1.4 (in principle 
only calibrated for N nn = 25). The dotted line corresponds to 
the number of clusters where substructure is detected along a 
given fraction of sightlines by the requirement P au t < 0.05. As 
the number of galaxies per test subgroup is changed from 10 or 
25, two other canonical values used in the literature, the fraction 
of sightlines per cluster with substructure decreases. The case 
shown here tends to have more substructure than these other 
two choices, because most of the clusters have the square root 
of number of galaxies < 10, and substructure detections seem to 
decrease as N n yi increases. 



structure criterion. Substructure detection using the require- 
ment of P su b < 0.05 occurred more often than detection 
using the requirement A/N nn > 1.4 (and almost always de- 
tected substructure when the latter test did). As the number 
of nearest neighbor galaxies used in the finder increased (i.e. 
larger substructures were sought), substructure was found 
along fewer lines of sight. Conversely, as the luminosity 
cut was lowered to include more galaxies, substructure was 
found more often. Lines of sight with substructure were not 
in 1-1 correspondence for different luminosity cuts, but the 
fraction of lines of sight with substructure was correlated 
cluster by cluster as the luminosity cut changed. 

The number of lines of sight with substructure per 
cluster was correlated with the number of galaxies in the 
largest substructure, and the number of galaxies in sub- 
groups with > 8 galaxies. (This is likely also related 
to the increas e of substructure found for recent merg- 
ers noted in lEspino-Briones. Plionis fc Ragone-Figueroa ' 



l|2007h : iRagone-Figueroa fc Plionisl \200H ) and mentioned 
earlier.) More substructure was found in higher mass clus- 
ters as well. To try to minimize the mass dependence, clus- 
ters were divided into 6 mass bins and substructure lines of 
sight were counted for both the top and bottom quartiles 
of galaxies in large subgroups. In these bins the clusters 
with more galaxies in large subgroups tended to have more 
lines of sight with substructure (but not always, and the 
results varied with substructure test criteria). If only sub- 
groups which had fallen in recently (i.e. within 0.6 Gyr) were 



6 SUMMARY 

This note considered subgroups in clusters, i.e. galaxies 
which shared the same halo before they fell into the galaxy 
cluster, in a relatively large volume and high resolution N- 
body simulation. The cluster sample had 243 clusters at 
2 = 0.1. Properties of the distribution of subgroup popu- 
lations, and how their orientations relative to line of sight of 
observation affected a variety of cluster mass measurements 
were considered, along with one test of cluster substructure. 

A significant fraction of cluster galaxies are in these sub- 
groups, with larger subgroups more likely in more massive 
(i.e. richer) clusters. In roughly half of the simulation clus- 
ters with M > 1O 14 /i _1 M0, at least one large (> 8 galaxies) 
subgroup is present, and 15% of the clusters have more than 
one large subgroup. The large subgroups tend to have higher 
galaxy densities and smaller rms velocities than the full sam- 
ple of galaxies in their host cluster. These correlations often 
remain many Gyr after infall. The largest subgroups tend to 
be found along the cluster long axis, with average velocities 
tending to be directed along this axis as well. Observation- 
ally, cluster mass measurements on average increase when 
the largest cluster subgroup lies along the line of sight, clus- 
ter by cluster, for two richness mass measurements and es- 
pecially for velocity dispersions. The larger the cluster sub- 
group, the more likely substructure will be detected by the 
Dressler-Shectman test, but the likelihood of detection also 
depends upon the line of sight direction, the number of near- 
est neighbors used in the test, and the luminosity cut for the 
galaxies included. 

It would be interesting to see how these subgroups 
evolve compared to the gas of their original infall halos, 
and to understand what other gala xy properties (besides 
merge rs and those dis cusse d i n iKnebe. Gill . 
' 20041'); iKnebe et all (120061 1; iMcIntosh ei"air 



Angulo et all J2009I); (Li fc Md (1200911: ISimha et aj 
Wetzel. Cohn fc Whitd (|2009h : iKlimentowski et all 



depend upon 
galaxies. 



subgroup membership for 
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